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Abstract Reliable predictions of general relativity theory are extracted using ap- 
proximation methods. Among these, the powerful post-Newtonian approximation 
provides us with our best insights into the problems of motion and gravitational 
radiation of systems of compact objects. This approximation has reached an im- 
pressive mature status, because of important progress regarding its theoretical foun- 
dations, and the successful construction of templates of gravitational waves emitted 
by inspiralling compact binaries. The post-Newtonian predictions are routinely used 
for searching and analyzing the very weak signals of gravitational waves in current 
generations of detectors. High-accuracy comparisons with the results of numerical 
simulations for the merger and ring-down of binary black holes are going on. In 
this article we give an overview on the general formulation of the post-Newtonian 
approximation and present up-to-date results for the templates of compact binary 
inspiral. 



1 Introduction 

Although relativists admire the mathematical coherence — and therefore beauty 
— of Einstein's general relativity, this theory is not easy to manage when draw- 
ing firm predictions for the outcome of laboratory experiments and astronomical 
observations. Indeed only few exact solutions of the Einstein field equations are 
known, and one is obliged in most cases to resort to approximation methods. It is 
beyond question that approximation methods in general relativity do work, and in 
some cases with some incredible precision. Many of the great successes of general 
relativity were in fact obtained using approximation methods. However because of 
the complexity of the field equations, such methods become awfully intricate at 
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high approximation orders. On the other hand it is difficult to set up a formalism 
in which the approximation method would be perfectly well-defined and based on 
clear premises. Sometimes it is impossible to relate the approximation method to 
the exact framework of the theory. In this case the only thing one can do is to rely 
on the approximation method as the only representation of real phenomena, and to 
discover "empirically" that the approximation works well. 

The most important approximation scheme in general relativity is the post- 
Newtonian expansion, which can be viewed as an expansion when the speed of 
light c tends to infinity, and is physically valid under the assumptions of weak gravi- 
tational field inside the source and of slow internal motion. The post-Newtonian ap- 
proximation makes sense only in the near-zone of an isolated matter source, defined 
as r <C A, where A = cT is the wavelength of the emitted gravitational radiation, 
with T being a characteristic time scale of variation of the source. The approxima- 
tion has been formalized in the early days of general relativity by Einstein 0"), De 
Sitter 12 13, Lorentz & Droste |4|. It was subsequently developed notably by Ein- 
stein, Infeld & Hoffmann [5,1, Fock EG), Plebanski & Bazanski Q, Chandrasekhar 
and collaborators JI] [TO] [H], Ehlers and his school 02] [13] O, Papapetrou and 
coworkers ifTBlfTBI . 

Several long-standing problems with the post-Newtonian approximation for gen- 
eral isolated slowly-moving systems have hindered progress untill recently. At high 
post-Newtonian ordners some divergent Poisson-type integrals appear, casting a 
doubt on the physical soundness of the approximation. Linked to that, the domain 
of validity of the post-Newtonian approximation is limited to the near-zone of the 
source, making it a priori difficult to incorporate into the scheme the condition of 
no-incoming radiation, to be imposed at past null infinity from an isolated source. 
In addition, from a mathematical point of view, we do not know what the "reliabil- 
ity" of the post-Newtonian series is, i.e. if it comes from the Taylor expansion of a 
family of exact solutions. 

The post-Newtonian approximation gives wonderful answers to the problems of 
motion and gravitational radiation, two of general relativity's corner stones. Three 
crucial applications are: 

1 . The motion of N point-like objects at the first post-Newtonian level J5] i.e. 1PN0 
is taken into account to describe the Solar System dynamics (motion of the cen- 
ters of mass of planets); 

2. The gravitational radiation-reaction force, which appears in the equations of mo- 
tion at the 2.5PN order lfT7l[T8l[T9ll20l . has been experimentally verified by the 
observation of the secular acceleration of the orbital motion of the Hulse-Taylor 
binary pulsar PSR 1913+16 ED EH E3; 

3. The analysis of gravitational waves emitted by inspiralling compact binaries — 
two neutron stars or black holes driven into coalescence by emission of gravi- 



1 As usual, we refer to «PN as the order equivalent to terms ~ (v/c) 2 " in the equations of motion be- 
yond the Newtonian acceleration, and in the asymptotic waveform beyond the Einstein quadrupole 
formula, where v denotes the binary's orbital velocity and c is the speed of light. 



Post-Newtonian theory and the two-body problem 



3 



tational radiation — necessitate the prior knowledge of the equations of motion 
and radiation field up to high post-Newtonian order. 

Strategies to detect and analyze the very weak signals from inspiralling compact 
binaries involve matched filtering of a set of accurate theoretical template wave- 
forms against the output of the detectors. Measurement-accuracy analyses have 
shown that, in order to get sufficiently accurate theoretical templates, one must at 
least include conservative post-Newtonian effects up to the 3PN level, and radiation- 
reaction effects up to 5.5PN order, i.e. 3PN beyond the leading order of radiation 
reaction which is 2.5PN ||24l |25] |26l |27l |28l |29l . 

The appropriate description of inspiralling compact binaries is by two structure- 
less point-particles, characterized solely by their masses m\ or 1112 (and possibly 
their spins). Indeed, most of the non-gravitational effects usually plaguing the dy- 
namics of binary star systems, such as the effects of a magnetic field, of an inter- 
stellar medium, the influence of the internal structure of the compact bodies, are 
dominated by purely gravitational effects. Inspiralling compact binaries are very 
clean systems which can essentially be described in pure general relativity. 

Although point-particles are ill-defined in the exact theory, they are admissible 
in post-Newtonian approximations. Furthermore the model of point particles can 
be pushed to high post-Newtonian order, where an a priori more realistic model 
involving the internal structure of compact bodies would fail through becoming in- 
tractable. However there is an important worry: a process of regularization to deal 
with the infinite self-field of point particles is crucially needed. The regularization 
should be carefully defined to be implemented at high orders. It should hopefully be 
followed by a renormalization. 

The orbit of inspiralling compact binaries can be considered to be circular, apart 
from the gradual inspiral, with an excellent approximation. Indeed, gravitational ra- 
diation reaction forces tend to circularize rapidly the orbital motion. At each instant 
during the gradual inspiral, the eccentricity e of the orbit is related to the instanta- 
neous frequency (0 = 2k/P by l30l 

e 2 ~ const co 19/9 (fore<l). (1) 

For instance one can check that the eccentricity of a system like the binary pulsar 
PSR 1913+16 will be e ~ 5 10~ 6 when the gravitational waves become visible by 
the detectors, i.e. when the signal frequency after its long chirp reaches /= co/n ~ 
10 Hz. Only systems formed very late, near their final coalescence, could have a 
non-negligible residual eccentricity. 

The intrinsic rotations or spins of the compact bodies could play an important 
role, yielding some relativistic spin-orbit and spin-spin couplings, both in the bi- 
nary's equations of motion and radiation field. The spin of a rotating body is of the 
order of S ~ ma v sp i n , where m and a denote the mass and typical size of the body, 
and where v sp i n represents the velocity of the body's surface. In the case of compact 
bodies we have a ~ Gm/c 2 , and for maximal rotation v sp i n ~ c; for such objects the 
magnitude of the spin is roughly S ~ Gm 2 /c. It is thus customary to introduce a 
dimensionless spin parameter, generally denoted by %, defined by 
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S 



(2) 



c 



We have % < 1 for black-holes, and % ^ 0.63 — 0.74 for neutron stars (depending 
on the equation of state of nuclear matter inside the neutron star). For binary pulsars 
such as PSR 1913+16, we have X < 5 10~ 3 . Considering models of evolution of ob- 
served binary pulsar systems when they become close to the coalescence we expect 
that the spins will make a negligible contribution to the accumulated phase in this 
case. However, astrophysical observations suggest that black holes can have non- 
negligible spins, due to spin up driven by accretion from a companion during some 
earlier phase of the binary evolution. For a few black holes surrounded by matter, 
observations indicate a significant intrinsic angular momentum and the spin may 
even be close to its maximal value. However, very little is known about the black- 
hole spin magnitudes in binary systems. For black holes rotating near the maximal 
value the templates of gravitational waves need to take into account the effects of 
spins, both for a successful detection and an accurate parameter estimation. 

We devote Part A of this article, including Sections [2]to [7] to a general overview 
of the formulation of the post-Newtonian approximation for isolated sources. We 
emphasize that the approximation can be carried out up to any post-Newtonian 
order, without the aforesaid problem of divergences. The main technique used is 
the matching of asymptotic expansions which permits to obtain a complete post- 
Newtonian solution incorporing the correct boundary conditions at infinity. Then 
Part B, i.e. Sections[8]to[T3] will deal with the application to systems composed of 
compact objects. The subtle issues linked with the self-field regularization of point- 
particles are discussed in some details. The results for the two-body equations of 
motion and radiation field at the state-of-the-art 3PN level are presented in the case 
of circular orbits appropriate to inspiralling compact binaries. We put the accent 
on the description of spinning particles and particularly on the spin-orbit coupling 
effect on the binary's internal energy and gravitational-wave flux. 



2 Einstein field equations 

General relativity is based on two independent tenets, the first one concerned with 
the dynamics of the gravitational field, the second one dealing with the coupling of 
all the matter fields with the gravitational field. Accordingly the action of general 
relativity is made of two terms, 
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gR + ^matter [^,8ap] ■ 



(3) 



The first term represents the kinetic Einstein-Hilbert action for gravity and tells that 
the gravitational field g a p propagates like a pure spin-2 field. Here R is the Ricci 
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scalar and g = det(g a p ) is the determinant of the metric. The second term expresses 
the fact that all matter fields (collectively denoted by f ) are minimally coupled to 
the metric g a p which defines the physical lengths and times as measured in local 
laboratory experiments. The field equations are obtained by varying the action with 
respect to the metric (such that 8g a p = when \xp | — > °°) and form a system of ten 
differential equations of second order, 

G a P[g,dg^g] = 8 -^T^[g}. (4) 

The Einstein tensor G a $ = R a $ — ^Rg a P is generated by the matter stress-energy 
tensor T a P = (2/ '\/— g)(5S ma tter/ ^gap)- F° ur equations give, via the contracted 
Bianchi identity, the conservation equation for the matter system as 

V M G^=0 =^> ^7^=0, (5) 

which must be solved conjointly with the Einstein field equations for the gravita- 
tional field. The matter equation (O reads also 

<?M (V=I#) = \^f—gT» v d agtlv . (6) 

Let us introduce an asymptotically Minkowskian coordinate system such that the 
gravitational-wave amplitude, defined by h a P = ^/—gg a P — rj^Q is divergence- 
less, i.e. satisfies the de Donder or harmonic gauge condition 

d^=0. (7) 

With this coordinate choice the Einstein field equations can be recast into the 
d' Alembertian equation 

Uh ^ = i^ T «£ . (8) 

c 4 

Here □ = rj^d^dy is the usual (/?af-spacetime) d' Alembertian operator. The source 
term T a P can rightly be interpreted as the "effective" stress-energy distribution of 
the matter and gravitational fields in harmonic coordinates; note that T a P is not a 
tensor and we shall often call it a pseudo-tensor. It is conserved in the sense that 

d^=0. (9) 

This is equivalent to the condition of harmonic coordinates (0 and to the covariant 
conservation §5$ of the matter tensor. The pseudo-tensor is made of the contribution 
of matter fields described by T a $ , and of the gravitational contribution A which 
is a complicated functional of the gravitational field variable h^ v and its first and 
second derivatives. Thus, 



2 Here, g a P denotes the contravariant metric, inverse of the covariant metric g a p, and r\ a P repre- 
sents an auxiliary Minkowski metric. We assume that our spatial coordinates are Cartesian so that 
fa"*) = fao0)=diag(-l, 1,1,1). 
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t a ^ = \g\T^ + -—A^[h,dh,d 2 h}. (10) 
lonG 

The point is that A a $ is at least quadratic in the field strength /j^ v , so the field 
equations (O are naturally amenable to a perturbative non-linear treatment. The 
general expression is 

-g^gvodph^d^P -g^g va d p h ao d,h v P +g^gP°d p h a »d a llP v 
+ l - {ig^g^ - g «V V ) (gpege* \goegp*) d^dyh" . (11) 

To select a physically sensible solution of the field equations in the case of a 
bounded matter system, we must impose a boundary condition at infinity, namely 
the famous no-incoming radiation condition, which ensures that the system is truly 
isolated from other bodies in the Universe. In principle the no-incoming radiation 
condition is to be formulated at past null infinity . Here, we can simplify the 
formulation by taking advantage of the presence of the Minkowski background r\ a $ 
to define the no-incoming radiation condition with respect to Minkowskian past null 
infinity say ^~ . Within approximate methods this is legitimate as we can view the 
gravitational field as propagating on the fiat background Tj a p ; indeed T] a p does exist 
at any finite order of approximation. 

The no-incoming radiation condition should be such that it suppresses at «/Jj~ any 
homogeneous (regular in R 4 ) solution of the d'Alembertian equation D/ihom = 
in a neighborhood of J?^ . We have at our disposal the Kirchhoff formula which 
expresses hh om at some field point (x',f') m terms of its values and its derivatives on 
a sphere centered on x' with radius p = |x' — x| and at retarded time t =t' — p/c, 



-T- (P^hom) + - -=r {phhom) 
dp c dt 



(x,f), (12) 



where dQ is the solid angle spanned by the unit direction N = (x — x')/p. From 
this formula we deduce the no-incoming radiation condition as the following limit 



at JF- , i.e. when r = |x| — >• +°° with t + r/c= constH 



lim 



dr c dt 



(13) 



Now if h a P satisfies ( fT3b . so does the pseudo-tensor built on it, and then it is 
clear that the retarded integral of T ^ does satisfy the same condition. Thus we infer 
that the unique solution of the Einstein equation (H)) reads 



3 In fact we obtain also the auxiliary condition that rdph a P should be bounded near J 2 ^ . This 
comes from the fact that p differs from r and we have p = r — x' ■ n+ ff(l/r) with n = x/r. 
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where the retarded integral takes the standard form 



(□ 



1 f dV 



x a P(x',t-\x-x'\/c). 



(15) 



'r 



An J |x-x'| 



Notice that since T a P depends on h^ v and its derivatives, the equation (fl4l i is to 
be viewed as an integro-differential equation equivalent to the Einstein equation I© 
with no-incoming radiation condition. 



3 Post-Newtonian iteration in the near zone 

In this Section we proceed with the post-Newtonian iteration of the field equations 
in harmonic coordinates in the near zone of an isolated matter distribution. We have 
in mind a general hydrodynamical fluid, whose stress-energy tensor is smooth, i.e. 
T a P G C°°(IR 4 ). Thus the scheme a priori excludes the presence of singularities; 
these will be dealt with in later Sections. 

Let us remind that the post-Newtonian approximation in "standard" form (e.g. 
BT] [13] [14)) is plagued with some apparently inherent difficulties, which crop up at 
some high post-Newtonian order like 3PN. Up to the 2.5PN order the approxima- 
tion can be worked out without problems, and at the 3PN order the problems can 
generally be solved for each case at hand; but the problems worsen at higher or- 
ders. Historically these difficulties, even appearing at higher approximations, have 
cast a doubt on the actual soundness, from a theoretical point of view, of the post- 
Newtonian expansion. Practically speaking, they posed the question of the reliabil- 
ity of the approximation, when comparing the theory's predictions with very precise 
experimental results. In this Section and the next one we assess the nature of these 
difficulties — are they purely technical or linked with some fundamental drawback 
of the approximation scheme? — and eventually resolve them. 

We first distinguish the problem of divergences in the post-Newtonian expan- 
sion: in higher approximations some divergent Poisson-type integrals appear. Re- 
call that the post-Newtonian expansion replaces the resolution of an hyperbolic- 
like d' Alembertian equation by a perturbatively equivalent hierarchy of elliptic-like 
Poisson equations. Rapidly it is found during the post-Newtonian iteration that the 
right-hand-sides of the Poisson equations acquire a non-compact support (it is dis- 
tributed all over space), and that the standard Poisson integral diverges because of 
the bound of the integral at spatial infinity, i.e. r = |x| — s- +°°, with t = const. 

The divergencies are linked to the fact that the post-Newtonian expansion is actu- 
ally a singular perturbation, in the sense that the coefficients of the successive pow- 
ers of 1 jc are not uniformly valid in space, typically blowing up at spatial infinity 
like some positive powers of r. We know for instance that the post-Newtonian ex- 
pansion cannot be "asymptotically flat" starting at the 2PN or 3PN level, depending 
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on the adopted coordinate system 1321 , The result is that the standard Poisson inte- 
grals are in general badly-behaving at infinity. Trying to solve the post-Newtonian 
equations by means of the Poisson integral does not a priori make sense. This does 
not mean that there are no solutions to the problem, but simply that the Poisson in- 
tegral does not constitute the good solution of the Poisson equation in the context 
of post-Newtonian expansions. So the difficulty is purely of a technical nature, and 
will be solved once we succeed in finding the appropriate solution to the Poisson 
equation. 

We shall now prove (following [33)) that the post-Newtonian expansion can be 
indefinitely iterated without divergences^ Let us denote by means of an overline 
the formal (infinite) post-Newtonian expansion of the field inside the source's near- 
zone, which is of the form 

h aP (x,t,c) = f t ^h a P(x,t,lnc). (16) 

^ C n 

n=2 

The «-th post-Newtonian coefficient is naturally the factor of the «-th power of 1 /c; 
however, we know (36| that the post-Newtonian expansion also involves some log- 
arithms of c, included for convenience here into the definition of the coefficients /;„. 
For the stress-energy pseudo-tensor (TToT > we have the same type of expansion, 

f^(x,/,c)= £ \z a ^(x,tMc), (17) 

n=—2 

The expansion starts with a term of order c 2 corresponding to the rest mass-energy 
(x a P has the dimension of an energy density). Here we shall always understand the 
infinite sums such as (fTo*ll-(fT7]i in the sense of formal series, i.e. merely as an ordered 
collection of coefficients. Because of our consideration of regular extended matter 
distributions the post-Newtonian coefficients are smooth functions of space-time. 

Inserting the post-Newtonian ansatz into the Einstein field equation ((HJ and 
equating together the powers of 1/c results is an infinite set of Poisson-type equa- 
tions (Vn > 2), 

Ah a ^ = l6nG T a P + d?h a P. (18) 

" n-4 n-2 

The second term comes from the split of the d'Alembertian into a Laplacian and a 
second time derivative: □ = A — This term is zero when n —2 and 3. We pro- 
ceed by induction, i.e. fix some post-Newtonian order n, assume that we succeeded 
in constructing the sequence of previous coefficients p h for p<n—l, and from this 
show how to infer the next-order coefficient „/?. 

To cure the problem of divergencies we introduce a generalized solution of the 
Poisson equation with non-compact support source, in the form of an appropriate 
finite part of the usual Poisson integral obtained by regularization of the bound at 



4 An alternative solution to the problem of divergencies, proposed in 1341 [35), is based on an 
initial-value formalism, which avoids the appearance of divergencies because of the finiteness of 
the integration region. 
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infinity by means of a specific process of analytic continuation. For any source term 
like „T, we multiply it by the "regularization" factor 



Ixf 



(19) 



where B £ C is a complex number and ro denotes an arbitrary length scale. Only then 
do we apply the Poisson integral, which therefore defines a certain function of B. 
The well-definedness of that integral heavily relies on the behavior of the integrand 
at the bound at infinity. There is no problem with the vicinity of the origin inside 
the source because of the smoothness of the pseudo-tensor. Then one can prove 1331 
that the latter function of B generates a (unique) analytic continuation down to a 
neighborhood of the value of interest B = 0, except at B = itself, at which value 
it admits a Laurent expansion with multiple poles up to some finite order. Then, we 
consider the Laurent expansion of that function when B — > and pick up the finite 
part, or coefficient of the zero-th power of B, of that expansion. This defines our 
generalized Poisson integral: 

A -1 [ T™' 5 ] (x,/) = — FP [ -^L\x'\ B T a P(x',t). (20) 

1 n 1 Alt B=Q J |x — X | n 

The integral extends over all three-dimensional space but with the latter finite-part 
regularization at infinity denoted FPg = o. The main properties of our generalized 
Poisson operator is that it does solve the Poisson equation, namely 



i [t 0! ^)=t^, (21) 



A (a~ 

and that the so defined solution A~ x „x owns the same properties as its source „T, 
i.e. the smoothness and the same type of behavior at infinity. 

The most general solution of the Poisson equation ( f]~8T > will be obtained by appli- 
cation of the previous generalized Poisson operator to the right-hand-side of ( [PHI ), 
and augmented by the most general homogeneous solution of the Poisson equation. 
Thus, we can write 

h a P = 16nGA- l [ l^]+d?A- l [Ti 2 a P]+Y,Bf(t)x L . (22) 



1=0' 



The last term represents the general solution of the Laplace equation which is regular 
at the origin r = |x| = 0. It can be written, using the symmetric-trace-free (STF) 
language, as a multipolar series of terms of the type i^jj and multiplied by some 



5 Here L = ;'i • • -i( denotes a multi-index composed of i multipolar spatial indices i\ , ■ ■ ■ ,ig (ranging 
from 1 to 3); xl = x- n ■ ■ ■xu is the product of i spatial vectors x,; di = • • -(3, f is the product of 
I partial derivatives d, = d / dx' ; in the case of summed-up (dummy) multi-indices L, we do not 
write the I summations from 1 to 3 over their indices; the STF projection is indicated with a hat, 
i.e. xi = STF[xJ and similarly ot = STF[ot], or sometimes using brackets surrounding the indices, 
x<l> =x L . 
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STF-tensorial functions of time „Bi(t). These functions will be associated with the 
radiation reaction of the field onto the source; they will depend on which boundary 
conditions are to be imposed on the gravitational field at infinity from the source. 

It is now trivial to iterate the process. We substitute for n -?h in the right-hand- 
side of (l22l the same expression but with n replaced by n — 2, and similarly come 
down until we stop at either one of the coefficients q!i = or \h = 0. At this point 
,,h is expressed in terms of the "previous" p t's and p Bls with p<n — 2. To finalize 
the process we introduce what we call the operator of the "instantaneous" potentials 
□ f 1 . Our notation is chosen to contrast with the standard operators of the retarded 
and advanced potentials D^ 1 and D^ 1 , see (fl31 l. However, beware of the fact that 
unlike A the operator D^ 1 will be defined only when acting on a post-Newtonian 
series such as T. Indeed, we pose 

where A~ k ~ l is the k-th iteration of the operator ( f20b . It is readily checked that in 
this way we have a solution of the source-free d' Alembertian equation, 

□ (□f 1 ^]) =T a K (24) 

On the other hand, the homogeneous solution in ( 1221 will yield by iteration an ho- 
mogeneous solution of the d' Alembertian equation which is necessarily regular at 
the origin. Hence it should be of the anti-symmetric type, i.e. be made of the dif- 
ference between a retarded multipolar wave and the corresponding advanced wave. 
We shall therefore introduce a new definition for some STF-tensorial functions A^(f ) 
parametrizing those advanced-minus-retarded free waves. It will not be difficult to 
relate the post-Newtonian expansion of Ai(t) to the functions „Bl(0 which were 
introduced in d22l . Finally the most general post-Newtonian solution, iterated ad 
infinitum and without any divergences, is obtained into the form 

-ap_\6%G lf p , 4G f A?(t-r/c)-A?(t + r/c) 

We shall refer to the Al(/)'s as the radiation-reaction functions. If we stay at the 
level of the post-Newtonian iteration which is confined into the near zone we can- 
not do more than (l25l i: there is no means to compute the radiation-reaction func- 
tions Ai(t). We are here touching the second problem faced by the standard post- 
Newtonian approximation. 
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4 Post-Newtonian expansion calculated by matching 

As we now understand this problem is that of the limitation to the near zone. Indeed 
the post-Newtonian expansion assumes that all retardations r/c are small, so it can 
be viewed as a formal near-zone expansion when r — > 0, valid only in the region 
surrounding the source that is of small extent with respect to the wavelength of 
the emitted radiation: r <^ X. As we have seen, a consequence is that the post- 
Newtonian coefficients blow up at infinity, when r +°°. It is thus not possible, 
a priori, to implement within the post-Newtonian scheme the physical information 
that the matter system is isolated from the rest of the Universe. The no-incoming 
radiation condition imposed at past null infinity cannot be taken into account, 
a priori, within the scheme. 

The near-zone limitation can be circumvented to the lowest post-Newtonian or- 
ders by considering retarded integrals that are formally expanded when c — > +°° as 
series of "instantaneous" Poisson-like integrals iFJTl . This procedure works well up 
to the 2.5PN level and has been shown to correctly fix the dominant radiation reac- 
tion term at the 2.5PN order lfl~3l IT41 . Unfortunately such a procedure assumes fun- 
damentally that the gravitational field, after expansion of all retardations r/c —> 0, 
depends on the state of the source at a single time t, in keeping with the instan- 
taneous character of the Newtonian interaction. However, we know that the post- 
Newtonian field (as well as the source's dynamics) will cease at some stage to be 
given by a functional of the source parameters at a single time, because of the im- 
print of gravitational-wave tails in the near zone field, in the form of some modifi- 
cation of the radiation reaction force at the 1.5PN relative order ll37l[38l . Since the 
reaction force is itself of order 2.5PN this means that the formal post-Newtonian 
expansion of retarded Green functions is no longer valid starting at the 4PN order. 

The solution of the problem resides in the matching of the near-zone field to the 
exterior field, a solution of the vacuum equations outside the source which has been 
developed in previous works using some post-Minkowskian and multipolar expan- 
sions |[36l [39l . In the case of post-Newtonian sources, the near zone, i.e. r<A, 
covers entirely the source, because the source's radius itself is such that a <C X. 
Thus the near zone overlaps with the exterior zone where the multipole expansion is 
valid. Matching together the post-Newtonian and multipolar-post-Minkowskian so- 
lutions in this overlapping region is an application of the method of matched asymp- 
totic expansions, which has frequently been applied in the present context, both for 
radiation-reaction j40]|4T][37l[2ll an d wave-generation (42l [43] 04) problems. 

In the previous Section we obtained the most general solution d25l ) for the post- 
Newtonian expansion, as parametrized by the set of unknown radiation-reaction 
functions A/, (f). We shall now impose the matching condition 

J((h afi ) = Ji({h a P), (26) 

telling that the multipole decomposition of the post-Newtonian expansion h of the 
inner field, agrees with the near-zone expansion of the multipole expansion j£(h) 
of the external field. Here the calligraphic letter stands for the multipole de- 
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composition or far-zone expansion, while the overbar denotes the post-Newtonian 
or near-zone expansion. The matching equation results from the numerical equality 
h = j% (h), clearly verified in the exterior part of the near-zone, namely our over- 
lapping region a < r <C A. The left-hand-side is expanded when r — > +°° yielding 

(h) while the right-hand-side is expanded when r->0 leading to ^(h). The 
matching equation is thus physically justified only for post-Newtonian sources, for 
which the exterior near-zone exists. It is actually afunctional identity; it identifies, 
term-by-term, two asymptotic expansions, each of them being formally taken out- 
side its own domain of validity. In the present context, the matching equation insists 
that the infinite far-zone expansion (r — > °o) of the inner post-Newtonian field is 
identical to the infinite near-zone expansion (r — > 0) of the exterior multipolar field. 
Let us now state that d26l i, plus the condition of no-incoming radiation, permits de- 
termining all the unknowns of the problem: i.e., at once, the external multipolar 
decomposition -#(/z) and the radiation-reaction functions Al and hence the inner 
post-Newtonian expansion h. 

When applied to a multipole expansion such as that of the pseudo-tensor, i.e. 
.<#(f a ^), we have to define a special type of generalized inverse d'Alembertian 
operator, built on the standard retarded integral ( fT5l ). viz 

0^[^(^)}(x,t)^-^wJ ^L\^^(^)(x',t-\x-^\/c). (27) 

Like (n~5b this integral extends over the whole three-dimensional space, but a regu- 
larization factor |x'| B given by (fT9b has been "artificially" introduced for application 
of the finite part operation FPb = o. The reason for introducing such regularization is 
to cure the divergencies of the integral when |x'| — > 0; these are coming from the 
fact that the multipolar expansion is singular at the origin. We notice that this reg- 
ularization factor is the same as the one entering the generalized Poisson integral 
d20l i, however its role is different, as it takes care of the bound at |x'| =0 rather than 
at infinity. We easily find that this new object is a particular retarded solution of the 
wave equation, 

□ (□r'[^(t^)]) = ..#(T a ^). (28) 

Therefore ^(h a P) should be given by that solution plus a retarded homogeneous 
solution of the d'Alembertian equation (imposing the no-incoming radiation condi- 
tion), i.e. be of the type 

^) = l^n R -.[^)] «£(z^j5a^/£)\. (29) 

c c (=Q i. ^ r j 

Now the matching equation d26l i will determine both the Al's in d25l l and the Fl's in 
d29b . We summarize the results which have been obtained in l33l . 

The functions £/,(/) will play an important role in the following, because they 
appear as the multipole moments of a general post-Newtonian source as seen from 
its exterior near zone. Their closed-form expression obtained by matching reads 
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F^(t) = wJ d\\x\ B x L J^dzdtiz)^ (x,t-z\x\/c) . (30) 

Again the integral extends over all space but the bound at infinity (where the post- 
Newtonian expansion becomes singular) is regularized by means of the same finite 
part. The z-integration involves a weighting function 8?(z) defined by 

The integral of that function is normalized to one: dzd((z) = 1. Furthermore it 
approaches the Dirac function in the limit of large multipoles: lim^,*, 8((z) = 8(z). 
The multipole moments (f30b are physically valid only for post-Newtonian sources. 
As such, they must be considered only in a perturbative post-Newtonian sense. With 
the result ( f30b the multipole expansion (f29t is fully determined and will be exploited 
in the next Section. 

Concerning the near-zone field d25b we find that the radiation-reaction functions 
Al are composed of the multipole moments Fl which will also characterize "linear- 
order" radiation reaction effects starting at 2.5PN order, and of an extra contribution 
Rl which will be due to non-linear effets in the radiation reaction and turn out to 
arise at 4PN order. Thus, 

A?=Fp+R¥, (32) 

where Fl is given by (f30b and where Rl is defined from the multipole expansion of 
the pseudo-tensor as 

Rf(t) = FP / d 3 x\x\ B x L [ + ~dzy e (z)^(T a h (x,t-z\x\/c) . (33) 

Here the regularization deals with the bound of the integral at |x| =0. Since the 
variable z extends up to infinity these functions truly depend on the whole past- 
history of the source. The weighting function therein is simply given by Ye(z) = 
—2d((z), this definition being motivated by the fact that the integral of that function 
is normalized to one: J^°°dzYe(z) = 1@ The specific contributions due to R L in the 
post-Newtonian metric d25l l are associated with tails of waves 13711381 . The fact that 
the external multipolar expansion (t) is the source term for the function Rl, and 
therefore will enter the expression of the near-zone metric ( 125) . is a result of the 
matching condition d26l i and reflects of course the no-incoming radiation condition 
imposed at J*^ . 

The post-Newtonian metric dZSb is now fully determined. However let us now 
derive an interesting alternative formulation of it 1451 . To this end we introduce still 
another object which will be made of the expansion of the standard retarded integral 
( fTBT ) when c — > °o, but acting on a post-Newtonian source term T, 



6 This integral is a priori divergent, however its value can be obtained by invoking complex analytic 
continuation in ( 6 C. 
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(x,r) ee -J- £ U! f-LY FP / dV |xf Ix-x'l- 1 T^(x' )f 



4tt ^ n! \cdt J b=o 

(34) 

Each of the terms is regularized by means of the finite part to deal with the bound 
at infinity where the post-Newtonian expansion is singular. This regularization is 
crucial and the object should carefully be distinguished from the "global" solution 
□k 1 [t] defined by ( fT4l and in which the pseudo-tensor is not expanded in post- 
Newtonian fashion. We emphasize that d34l ) constitutes merely the definition of a 
(formal) post-Newtonian expansion, each term of which being built from the post- 
Newtonian expansion of the pseudo-tensor. Such a definition is of interest because it 
corresponds to what one would intuitively think as the "natural" way of performing 
the post-Newtonian iteration, i.e. by Taylor expanding the retardations as in iTJTIl . 
Moreover, each of the terms of the series d34l ) is mathematically well-defined thanks 
to the finite part, and can therefore be implemented in practical computations. The 
point is that d34l > solves, in a post-Newtonian sense, the wave equation, 

□ (□^pr " 5 ]) =t"^, (35) 



so constitutes a good prescription for a particular solution of the wave equation 
— as legitimate a prescription as ( f23l . Therefore d23l and ( 1341 should differ by 
an homogeneous solution of the wave equation which is necessarily of the anti- 
symmetric type. Detailed investigations yield 

- T'^i -^i^^{ f °'"- r/c) 2 r"" +r/c) } • 

(36) 

in which the homogeneous solution is parametrized precisely by the multipole- 
moment functions Fi{t). This formula is the basis of our writing of the new form of 
the post-Newtonian expansion. Indeed, by combining (l25l l and d36l l. we nicely get 

1* - ig° □;■ [r*] - g jffi { Rf( ' ~ " C) - g "' ( ' + " C) } ■ (37) 

which is our final expression for the general solution of the post-Newtonian field 
in the near-zone of any isolated matter distribution. This expression is probably the 
most convenient and fruitful when doing practical applications. 

We recognize in the first term of (l37T i (notwithstanding the finite part therein) 
the old way of performing the post-Newtonian expansion as it was advocated by 
Anderson & DeCanio OTI . For computations limited to the 3.5PN order, i.e. up 
to the level of the 1PN correction to the radiation reaction force, such first term 
is sufficient. However, at the 4PN order there is a fundamental breakdown of this 
scheme and it becomes necessary to take into account the second term in < [37T > which 
corresponds to non-linear radiation reaction effects associated with tails. 



Post-Newtonian theory and the two-body problem 



15 



Note that the post-Newtonian solution, in either form d25l l or (137] !, has been ob- 
tained without imposing the condition of harmonic coordinates in an explicit way, 
see (0. We have simply matched together the post-Newtonian and multipolar ex- 
pansions, satisfying the "relaxed" Einstein field equations ([8]) in their respective do- 
mains, and found that the matching determines uniquely the solution. An important 
check (carried out in ll33l l45ll ) is therefore to verify that the harmonic coordinate 
condition (0 is indeed satisfied as a consequence of the conservation of the pseudo- 
tensor (O, so that we really grasp a solution of the full Einstein field equations. 



5 Multipole moments of a post-Newtonian source 

The multipole expansion of the field outside a general post-Newtonian source has 
been obtained in the previous Section as[3 

where the multipole moments are explicitly given by (f30b . and the second piece 
reflects the non-linearities of the Einstein field equations and reads 

u ttP = n^[Jt(A a l > )]. (39) 

To write the latter expression we have used the fact that since the matter tensor T a P 
has a spatially compact support we have ^£(T a P ) = 0. Thus u a P is indeed generated 
by the non-linear gravitational source term ( fTTT i. We notice that the divergence of this 
piece, say w a = d^u 01 ^, is a retarded homogeneous solution of the wave equation, 
i.e. of the same type as the first term in d38l ). Now from w a we can construct a 
secondary object v a P which is also a retarded homogeneous solution of the wave 
equation, and furthermore whose divergence satisfies d^v"^ = — w a , so that it will 
cancel the divergence of u a P (see l44l for details). With the above construction of 
v a fi we are able to define the following combination, 

which will constitute the linearized approximation to the multipolar expansion 
^((h a $ ) outside the source. Then we have 

Jt(h a V) = Gti& + u a P + v a P . (41) 



7 An alternative formulation of the multipole expansion for a post-Newtonian source, with non- 
STF multipole moments, has been developed by Will and collaborators I46ll47ll48l . 



16 Luc Blanchet 

Having singled out such linearized part, it is clear that the sum of the second and 
third terms should represent the non-linearities in the external field. If we index 
those non-linearities by Newton's constant G, then we can prove indeed that u a P + 
v a P = ff(G 2 ). More precisely we can decompose u a P + v a P as a complete non- 
linearity or "post-Minkowskian" expansion of the type 

M ^+v^=£ CTh$y (42) 

m=2 

One can effectively define a post-Minkowskian "algorithm" l36l l44l able to con- 
struct the non-linear series up to any post-Minkowskian order m. The post-Minkows- 
kian expansion represents the most general solution of the Einstein field equations 
in harmonic coordinates valid in the vacuum region outside an isolated source. 

The above linearized approximation hry solves the linearized vacuum Einstein 
field equations in harmonic coordinates and it can be decomposed into multipole 
moments in a standard way |49l . Modulo an infinitesimal gauge transformation pre- 
serving the harmonic gauge, namely 

= + d a tf 1} + dPtpfa 77 a/3 <?M<o > (43) 
where the infinitesimal gauge vector pSL satisfies = 0, we can decompos^l 

•Si = + <ti*^> ' ,44b) 

% - - ? e (iigu) + } ■ <«<=> 

This decomposition defines two types of multipole moments, both assumed to be 
STF: the mass-type Il(w) and the current-type J/,(k). These moments can be arbi- 
trary functions of the retarded time u = t — r/c, except that the monopole and dipole 
moments (having I < 1) satisfy standard conservation laws, namely 

T (1) =I (2) =J (1) =0 . (45) 

The gauge transformation vector admits a decomposion in similar fashion, 



t ^ t. 



(46a) 



8 The supersript (k) refers to k time derivatives of the moments; e„/, t . is the Levi-Civita antisymmet- 
ric symbol such that £123 = 1. From here on the spatial indices such as j, ... will be raised and 
lowered with the Kronecker metric &/. They will be located lower or upper depending on context. 
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?EtK' (; y -0 + ZTT e -^- 1 (V 1 )} ' (46c) 



The six sets of STF multipole moments II, Jl, Wl, Xl, Yl and Zz, will collectively 
be called the multipole moments of the source. They contain the full physical in- 
formation about any isolated source as seen from its exterior near zone. Actually 
it should be clear that the main moments are li and Jl because the other moments 
Wl, • • • , Zi parametrize a linear gauge transformation and thus have no physical 
implications at the linearized order. However because the theory is covariant with 
respect to non-linear diffeomorphisms and not merely with respect to linear gauge 
transformations, the moments Wl, • • • , Zl do play a physical role starting at the 
non-linear level. We shall occasionaly refer to the moments Wl, Xl, Yl and Zl as 
the gauge moments. 

To express in the best way the source multipole moments, we introduce the fol- 
lowing notation for combinations of components of the pseudo-tensor f a P , 



(47a) 



- T°< 

Ei= — , (47b) 

c 

Z,j=T ij , (47c) 

where f" = 5,/f 7 . Here the overbar reminds us that we are exclusively dealing with 
post-Newtonian-expanded expressions, i.e. formal series of the type (1171 1. Then the 
general expressions of the "main" source multipole moments II and in the case 
of the time-varying moments for which £>2, are 

8 e+2 x ijL z\f\ , (48a) 



h(u) = FPy d 3 x\x\ B j dzEab^ 8 e x L _ l)a Z b 



c A {£+l){t + 2){2t + 5) 



c 2 {£ + 2)(2£ + 3) 



8 e+l x L _ l)ac E^y (48b) 



The integrands are computed at the spatial point x and at time m + z|x|/c, where u = 
t — r/cis the retarded time at which are evaluated the moments. We recall that z is the 
argument of the function 8( defined in OTb . Similarly we can write the expressions 
of the gauge-type moments W^, • • ■ , Z^. Notice that the source multipole moments 
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d48l ) have no invariant meaning; they are defined for the harmonic coordinate system 
we have chosen. 

Of what use are these results for the multipole moments 1^ and Jj,? From d44b 
these moments parametrize the linearized metric hi^, which is the "seed" of an in- 
finite post-Minkowskian algorithm symbolized by d42l i. For a specific application, 
i.e. for a specific choice of matter tensor like the one we shall describe in Section 
[H] the expressions d48b have to be worked out up to a given post-Newtonian order. 
The moments should then be inserted into the post-Minkowskian series (02) for the 
computation of the non-linearities. The result will be in the form of a non-linear 
multipole decomposition depending on the source moments 1^, J^, • • • , Z^, say 

JK{h a fi)= £ G m h«f n) {I L ,J L ,---}. (49) 

m=\ 

In the next Section we shall expand this metric at (retarded) infinity from the source 
in order to obtain the observables of the gravitational radiation field. 



6 Radiation field and polarisation waveforms 

The asymptotic waveform at future null infinity from an isolated source is the 
trans verse-traceless (TT) projection of the metric deviation at the leading order l/R 
in the distance R= |X| to the source, in a radiative coordinate system X^- = (cr,X)@ 
The waveform can be uniquely decomposed ll49l into radiative multipole compo- 
nents parametrized by mass-type moments Ul and current-type ones V/.. We shall 
define the radiative moments in such a way thay they agree with the f-th time deriva- 
tives of the source moments li and Jl at the linear level, i.e. 

U L = if + 0{G) , (50a) 
V L = jf + 0{G) . (50b) 

At the non-linear level the radiative moments will crucially differ from the source 
moments; the relations between these two types of moments will be discussed in the 
next Section. The radiative moments Ul(£/) and Vl(U) are functions of the retarded 
time U = T — R/c'm radiative coordinates. 

The asymptotic waveform at distance R and retarded time U is then given by 

AG + °° 1 f 2£ I 

h l T = ^ L \ N ^-2 UklL-2 ~ NaL-2 £ ab(k V l)bL - 2 } • (51) 



9 Radiative coordinates T and X, also called Bondi-type coordinates [501, are such that the met- 
ric coefficients admit an expansion when R — > +°° with U = T — R/c being constant, in simple 
powers of l/R, without the logarithms of R plaguing the harmonic coordinate system. Here U is 
a null or asymptotically null characteristic. It is known that the "far-zone" logarithms in harmonic 
coordinates can be removed order-by-order by going to radiative coordinates 1511 . 
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We denote by N = X/R = (Nj) the unit vector pointing from the source to the far- 
away detector. The TT projection operator reads S^ijki — &ik & 'ji ~ h &ki where 
= 8ij — NtNj is the projector orthogonal to the unit direction N. We project out 
the asymptotic waveform ( fBTI ) on polarization directions in a standard way. We de- 
note the two unit polarisation vectors by P and Q, which are orthogonal and trans- 
verse to the direction of propagation N (hence = P{Pj + QiQj). Our conventions 
and choice for P and Q will be specified in Section [T2l Then the two "plus" and 
"cross" polarisation states of the waveform are 



, _ PjPj QiQj , TT 

"H 2 ij ' 

_ P,Qj + PjQ, TT 

h x - 2 htj . 



(52a) 
(52b) 



Although the multipole decomposition dSTb entirely describes the waveform, 
it is also important, especially having in mind the comparison between the post- 
Newtonian results and numerical relativity 1521 . to consider separately the modes 
(£,m) of the waveform as defined with respect to a basis of spin-weighted spherical 
harmonics. To this end we decompose h + and h x as (see e.g. 15211531 ) 



(53) 



e=2m=-e 



where the spin-weighted spherical harmonics of weight —2 is a function of the 
spherical angles (©, <P) defining the direction of propagation N and reads 



'-2 



2£+l 

47T 



d fm (0)e 



i m <P 



k 2 



k=ki 



(ni 







cos ■ 



2(+m-2k-2 







sin ■ 



2k-m+2 



D 



lm _ (-)* ^(£ + m)\(£-m)\(£ + 2)l(l-2)l 
k\ {k-m + 2)\{l + m-k)\{l-k-2)\ 



(54a) 
(54b) 

(54c) 



HereA:i =max(0,m — 2) and&2 = min(£ + m,f — 2). Using the orthonormality prop- 
erties of these harmonics we obtain the separate modes hr m from the surface integral 
(with the overline denoting the complex conjugate) 



AO. h+-ih x Y im 2 (0,<f>). 



(55) 



On the other hand, we can also write //'" directly in terms of the radiative multipole 
moments Uz, and Vl, with result 



/; 



lm 



V2Rc e + 2 



u 



fin 



(56) 
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where U ft " and V*" are the radiative mass and current moments in non-STF guise. 
These are given in terms of the STF moments by 




£!V2(£ + l)(^lj a '" Vl ' (5?b) 

Here OC^ denotes the STF tensor connecting together the usual basis of spherical 
harmonics Y Cm to the set of STF tensors N L = STF(A^), recalling that Y lm and N L 
represent two basis of an irreducible representation of weight I of the rotation group. 
They are related by 



# L (0,<P)= £ ai m Y em (0,0), (58a) 

m=—i 

Y em (0,<P)= ( ^- a^N L (0,0), (58b) 



with the STF tensorial coefficient being 

af* = J d£2 N L Y bn . (59) 

The decomposition in spherical harmonic modes is especially useful if some of the 
radiative moments are known to higher post-Newtonian order than others. In this 
case the comparison with the numerical calculation ll52l l53l can be made for these 
individual modes with higher post-Newtonian accuracy. 



7 Radiative moments versus source moments 

The basis of our computation is the general solution of the Einstein field equa- 
tions outside an isolated matter system computed iteratively in the form of a post- 
Minkowskian or non-linearity expansion d49l (see details in ll36l[39l ). Here we give 
some results concerning the relation between the set of radiative moments {Ul, V/,} 
and the sets of source moments {II, Jl} and gauge moments {W^, • • • ,Zi}. Com- 
plete results up to 3PN order are available and have recently been used to control 
the 3PN waveform of compact binaries 1541 . 

Armed with definitions for all those moments, we proceed in a modular way. 
We express the radiative moments {Ul, V/.} in terms of some convenient intermedi- 
ate constructs {Ml,Sl} called the canonical moments. Essentially these canonical 
moments take into account the effect of the gauge transformation present in (l43l . 
Therefore they differ from the source moments {II, Jl} only at non-linear order. We 
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shall see that in terms of a post-Newtonian expansion the canonical and source mo- 
ments agree with each other up to 2PN order. The canonical moments are then con- 
nected to the actual source multipole moments {II, J/,} and {W^, • • • ,Z^}. The point 
of the above strategy is that the source moments (including gauge moments) admit 
closed-form expressions as integrals over the stress-energy distribution of matter 
and gravitational fields in the source, as shown in (|48| >. 

The mass quadrupole moment Uy (having i = 2) is known up to the 3PN order 
ll55ll . At that order it is made of quadratic and cubic non-linearities, and we have 



Uy(f/)=M|f(C/) + 
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(60) 



Notice the quadratic tail integral at 1.5PN order, the cubic tail-of-tail integral at 3PN 
order, and the non-linear memory integral at 2.5PN order 1 56 .57 .58 ,39 1 . The tail is 
composed of the coupling between the mass quadrupole moment Mjj and the mass 
monopole moment or total mass M; the tail-of-tail is a coupling between My and 
two monopoles M x M; the non-linear memory is a coupling My x M^. All these 
"hereditary" integrals imply a dependence of the waveform on the complete history 
of the source, from infinite past up to the current retarded time U = T — R/c. The 
constant uq in the tail integrals is defined by uq = ro/c, where ro is the arbitrary 
length scale introduced in < TT~9b . 

Note that the dominant hereditary integral is the tail arising at 1.5PN order in all 
radiative moments. For general £ we have at that order 
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where the constants Kg and Tig are given by 

2^ 2 + 5£ + 4 t?l 



(61a) 
(61b) 

(62a) 
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Now it can be proved that the retarded time U in radiative coordinates reads 

r 2GM ( r \ „ ( 1 
< V r 



U = t - («) 



where (t,r) are the harmonic coordinates. Inserting U into doTt we obtain the radia- 
tive moments expressed in terms of local source-rooted coordinates (t,r), e.g. 



. -u — r/c 



(64) 

This no longer depends on the constant uq — i.e. the «o gets replaced by the retar- 
dation time r/c. More generally it can be checked that mo always disappears from 
physical results at the end. On the other hand we can be convinced that the constant 
Kg (and Tig as well) depends on the choice of source -rooted coordinates (t,r). For in- 
stance if we change the harmonic coordinate system (t,r) to some "Schwarzschild- 
like" coordinates (t\r') such that t' = t and / = r + GM/c 2 , we get a new constant 
id = Ki + 1/2. Thus we have = 11/12 in harmonic coordinates [as shown in 
d60b l. but = 17/12 in Schwarzschild coordinates. 

We still have to relate the canonical moments {Ml,Sl} to the source multipole 
moments. As we said the difference between these two types of moments comes 
from the gauge transformation ( l46b and arises only at the small 2.5PN order. The 
consequence is that we have to worry about this difference only for high post- 
Newtonian waveforms. For the mass quadrupole moment My, the requisite cor- 
rection is given by 



W^Iy-WWlWj+^flV (65) 



4G 

M ij = l ij + - f 

where Iy denotes the source mass quadrupole, and where W is the monopole corre- 
sponding to the gauge moments (i.e. the moment having I — 0). Up to 3PN order, 
W is only needed at Newtonian order. The expression (f65t is valid in a mass-centred 
frame defined by the vanishing of the conserved mass dipole moment: I, = 0. 

Note that closed-form formulas generalizing (l65l l and similar expressions to all 
post-Newtonian orders (and all multipole interactions), if they exist are not known; 
they need to be investigated anew for specific cases. Thus it is convenient in the 
present approach to systematically keep the source moments {I/,, Jl,Wl,Xl, Yl,Zl} 
as the fundamental variables describing the source. 

B Inspiralling compact binaries 



8 Stress-energy tensor of spinning particles 



So far the post-Newtonian formalism has been developed for arbitrary matter dis- 
tributions. We want now to apply it to material systems made of compact objects 
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(neutron stars or black holes) which can be described with great precision by point 
masses. We thus discuss the modelling of point-particles possibly carrying some 
intrinsic rotation or spin. This means finding the appropriate stress-energy tensor 
which will have to be inserted into the general post-Newtonian formulas such as the 
expressions of the source moments ( |48T >. 

In the general case the stress-energy tensor will be the sum of a "monopo- 
lar" piece, which is a linear combination of monopole sources, i.e. made of Dirac 
delta-functions, plus the "dipolar" or spin piece, made of gradients of Dirac delta- 
functions. Hence we write 

T ali =TlL + T s %. (66) 

The monopole part takes the form of the stress-energy tensor for N particles (la- 
belled by A — 1 , • • • ,A0 without spin, reading in a four-dimensional picture 

T ap v [ + °° A „ (a p) 8^{x-y A ) 

Tmono = c L A Pa 11 a — / --, \ ■ (67) 

Here 5' 4 ' is the four-dimensional Dirac function. The world-line of particle A, de- 
noted y", is parametrized by the particle's proper time %a- The four-velocity is given 
by u" = dy"/dtA and is normalized to (^^ v )a m a m a = — c2 ' wnere (§hv)a denotes 
the metric at the particle's location. The four-vector p" is the particle's linear mo- 
mentum. For particles without spin we shall simply have p" = iiiau". However for 
spinning particles p" will differ from that and include some contributions from the 
spins as given by d75l l below. 

The dipolar part of the stress-energy tensor depends specifically on the spins, 
and reads in the classic formalism of spinning particles (due to Tulczyjew |59ll60ll . 
Trautman ED, Dixon E3, Bailey & Israel l63l ). 



•'spin — 2- V M 



A=l 



'dT^V^-^ 



(68) 



where V^, is the co variant derivative, and the anti-symmetric tensor S A represents 
the spin angular momentum of particle A. In this formalism the momentum-like 
quantity p" [entering d67l il is a time-like solution of the equation 

^ = (p>^<), (69) 

where D/dT^ denotes the covariant proper time derivative. The equation of trans- 
lational motion of the spinning particle, equivalent to the covariant conservation 
V^r afl = of the total stress-energy tensor d66l ). is the Papapetrou equation [64] 
involving a coupling to curvature, 

= 4«(*V)a- (70) 
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The Riemann tensor is evaluated at the particle's position A. This equation can also 
be derived directly from an action principle |63l . 

It is well-known that a choice must be made for a supplementary spin condition 
in order to fix unphysical degrees of freedom associated with an arbitrariness in the 
definition of the spin tensor S a P . This arbitrariness can be interpreted, in the case of 
extended bodies, as a freedom in the choice for the location of the center-of-mass 
worldline of the body, with respect to which the angular momentum is defined (see 
e.g. l65l ). An elegant spin condition is the covariant one, 

S a /p^ = 0, (71) 

which allows a natural definition of a spin four-(co)vector S„ such that 

sf = — _L= £ «/3mv 3_ s a (72) 

V-(S)A m AC 

For the spin vector itself, we can choose a four-vector which is purely spatial 
in the particle's instantaneous rest frame, where m" = (1,0). Therefore, we deduce 
that in any frame 

S£< = 0. (73) 

As a consequence of the covariant spin condition ( TTTb . we easily verify that the spin 
scalar is conserved along the trajectories, i.e. 

Sfs£ v = const. (74) 

Furthermore, we can check, using d7Tl i and also the law of motion d70l ). that the mass 
defined by m\c 2 = —p^P^ is also constant along the trajectories: tiia = const. Fi- 
nally, the relation linking the four-momentum and the four-velocity m" is readily 
deduced from the contraction of d69l ) with the four-momentum, which results in 

pi (pu) A + m 2 A c 2 u a A = l - S a / S v / u° (R^ avp ) A , (75) 

where (pu)^ = P^u^. Contracting further this relation with the four- velocity one 
deduces the expression of (pu)A and inserting this back into d75l ) yields the desired 
relation between p" and m". 

Focusing our attention on spin-orbit interactions, which are linear in the spins, 
we can neglect quadratic and higher spin corrections denoted £>'(S 2 )\ drastic sim- 
plifications of the formalism occur in this case. Since the right-hand-side of (l75l l is 
quadratic in the spins, we find that the four-momentum is linked to the four- velocity 
by the simple proportionality relation 

P a A=m A u a A + 0(f). (76) 



Hence, the spin condition ( TTTb becomes 
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S^uf l = ff{S 3 ). (77) 

Also, the equation of evolution for the spin, sometimes called the precessional equa- 
tion, follows immediately from the relationship (|69l together with the law (|76| > as 

Y ^ = ff{S 2 ) ^ = 0<s>). (78) 

dT A dT A 

Hence the spin vector S„ satisfies the equation of parallel transport, which means 
that it remains constant in a freely falling frame, as could have been expected be- 
forehand. Of course the norm of the spin vector is preserved along trajectories, 

S*S%= const. (79) 

In SectionQ~3]we shall apply this formalism to the study of spin-orbit effects in the 
equations of motion and energy flux of compact binaries. 



9 Hadamard regularization 

The stress-energy tensor of point masses has been defined in the previous Section by 
means of Dirac functions, and involves metric coefficients evaluated at the locations 
of the point particles, namely (gap) a- However, it is clear that the metric g a p be- 
comes singular at the particles. Indeed this is already true at Newtonian order where 
the potential generated by N particles reads 

w-i Gmn 

£/=£—, (80) 
s=i r B 

where r# = |x — yg\ is the distance between the field point and the particle B. Thus 
the values of U and hence of the metric coefficients [recall that goo = — 1 + 2U /c 2 + 
0{c ~ 4 )], are ill-defined at the locations of the particles. What we need is a self-field 
regularization, i.e. a prescription for removing the infinite self-field of the point 
masses. Arguably the choice of a particular regularization constitutes a fully quali- 
fied element of our physical modelling of compact objects. At Newtonian order the 
regularization of the potential (f80b should give the well-known result 

(C/).= I^, (81) 

B+A r ^B 

where r^B = \ya — Yb\ and the infinite self interaction term has simply been dis- 
carded from the summation. At high post-Newtonian orders the problem is not triv- 
ial and the self-field regularization must be properly defined. 

The post-Newtonian formalism reviewed in Sections[2H7]assumed from the start 
a continuous (smooth) matter distribution. Actually this formalism will be applica- 
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ble to singular point-mass sources, described by the stress-energy tensor of Section 
[8] provided that we suplement the scheme by the self-field regularization. Note that 
this regularization has nothing to do with the finite-part process FPg = o extensively 
used in the case of extended matter sources. The latter finite part was an ingredient 
of the rigorous derivation of the general post-Newtonian solution [see d20lil. while 
the self-field regularization is an assumption regarding a particular type of singular 
source. 

Our aim is to compute up to 3PN order the metric coefficients at the location 
of one of the particles: (g a p) A - At this stage different self-field regularizations are 
possible. We first review Hadamard's regularization [66, 67], which has proved to 
be very efficient for doing practical computations, but suffers from the important 
drawback of yielding some "ambiguity parameters", which cannot be determined 
within the regularization, starting at the 3PN order. 

Iterating the Einstein field equations with point-like matter sources (delta func- 
tions with spatial supports localized on Ja) yields a generic form of functions rep- 
resenting the metric coefficients in successive post-Newtonian approximations. The 
generic functions, say F(x), are smooth except at the points y^, around which they 
admit singular Laurent expansions in powers and inverse powers of r& = |x — | . 
When r^-s-Owe have (say, for any P <G N) 

fW=I^/»+^)- (82) 

p=Po A 

The coefficients a/p of the various powers of depend on the unit direction = 
(x — yx) jr A of approach to the singular point A. The powers p are relative integers, 
and are bounded from below by po <G Z. The Landau o-symbol for remainders takes 
its standard meaning. The a/p's depend also on the (coordinate) time ?, through 
their dependence on velocities Vfl(f) and relative positions ysc( f ) = ys(f) — yc( f ); 
however the time t is purely "spectator" in the regularization process, and thus will 
not be indicated. The coefficients a/p for which p < are referred to as the singular 
coefficients of F around A. 

The function F being given that way, we define the Hadamard partiefinie as the 
following value of F at the location of the particle 

(F) A = {f )= /^/ (n A ), (83) 

A J 47T A 

where &Q.A denotes the solid angle element centred on y A and sustained by n^. The 
brackets () mean the angular average. The second notion of Hadamard partie finie 
concerns the integral J d 3 xF, which is generically divergent at the points y A . Its 
partie finie (in short Pf) is defined by 

H S1 ... SN [d 3 xF = lim ( / d 3 xF+47r Y D A (s)\ . (84) 

J {M3\(JB A ( S ) a~1 J 



With this definition it is immediate to check that the previous Newtonian result jS 1 1 will hold. 
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The first term integrates over the domain M? \ Ua=i Ba( s ) defined as R 3 deprived 
from the N spherical balls Ba(s) = {x; r& < s} of radius s and centred on the points 
y^. The second term is the opposite of the sum of divergent parts associated with 
the first term around each of the particles in the limit where s->0. We have 



Da(s) = E 



4 s p+3 



(f P ) 



p + 3 A 



I - ) </-3) 
k S AJ a 



(85) 



Since the divergent parts are cancelled (by definition) the Hadamard partie finie is 
obtained in the limit s — > 0. Notice that as indicated in d84l ) the Hadamard partie-finie 
integral is not fully specified: it depends on N strictly positive and a priori arbitrary 
constants s\,--- ,syv parametrizing the logarithms in ( l85l ). 

We have seen that the post-Newtonian scheme consists of breaking the hyper- 
bolic d'Alembertian operator □ into the elliptic Laplacian A and the retardation 
term c~ 2 <9, 2 considered to be small, and put in the right-hand-side of the equation 
where it can be iterated; see (1181 . We thus have to deal with the regularization of 
Poisson integrals, or iterated Poisson integrals, of some generic function F . The 
Poisson integral will be divergent and we apply the prescription (184) , Thus, 



P(x') = -^Pf„ 



d 3 x 



(86) 



This definition is valid for each field point x' different from the Va's, and we want 
to investigate the singular limit when x' tends to one of the source points y^, so 
as to define the object (P)a- The definition ( T83T > is not directly applicable because 
the expansion of the Poisson integral P(x') when x' — » y^ will involve besides the 
normal powers of = |x' — va| some logarithms of r^. The proper way to define 
the Hadamard partie finie in this case is to include the lnr^ into the definition ([83} 
as if it were a mere constant parameter. With this definition we arrive at 



1 



(p), = -^pf„y-F( X ) 



d 3 x 



In ( ^- | -1 

sa 



(f-2). 
A 



(87) 



The first term is given by a partie-finie integral following the definition ([84); the 
second involves the logarithm of r* A . The constants s \ , ■ ■ ■ ,sn come from d86t . Since 
is actually tending to zero, In r* A represents a formally infinite "constant", which 
will ultimately parametrize the final Hadamard regularized 3PN equations of mo- 
tion. In the two-body case we shall find that the constants r* A are unphysical in the 
sense that they can be removed by a coordinate transformation [69|. Note that the 
apparent dependence of d87l l on the constant sa is illusory. Indeed the dependence 
on sa cancels out between the first and the second terms in the right-hand-side of 
([87), so the result depends only on r' A and the sg's for B ^ A. We thus have a simpler 



1 1 We consider only the local divergencies due to the singular points ■ The problem of divergen- 
cies of Poisson integrals at infinity is part of the general post-Newtonian formalism and has been 
treated in SectionfJ] 
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rewriting of (f87T > as 

( P )A = -^^- W /^W-</- a >- W 

Unfortunately, the constants sb for Z? 7^ A remaining in the result d88l will be the 
source of a genuine ambiguity. This ambiguity can in fact be traced back to the so- 
called "non-distributivity" of the Hadamard partie finie, a consequence of the pres- 
ence of the angular integration in (l83"T l. and implying that (FG)a 7^ (F)a (G)a in gen- 
eral. The non-distributivity arises precisely at the 3PN order both in the equations of 
motion and radiation field of point-mass binaries. At that order we are loosing with 
Hadamard's regularization an elementary rule of ordinary calculus. Consequently 
we expect that some basic symmetries of general relativity such as diffeomorphism 
invariance will be lost. However Hadamard's regularization can still be efficiently 
used to compute most of the terms in the equations of motion and radiation field 
at 3PN order; only a few ambiguous terms will show up which have then to be 
determined by another method. 



10 Dimensional regularization 

Dimensional regularization is an extremely powerful regularization which is free 
of ambiguities (at least up to the 3PN order). The main reason is that it is able to 
preserve the symmetries of classical general relativity; in fact dimensional regular- 
ization was invented ||701|7T1 as a means to preserve the gauge symmetry of pertur- 
bative quantum field theories. In the present context we shall show that dimensional 
regularization permits to resolve the problem of ambiguities arising at the 3PN or- 
der in Hadamard's regularization. We shall employ dimensional regularization not 
merely as a trick to compute some particular integrals which would otherwise be di- 
vergent or ambiguous, but as a fundamental tool for solving in a consistent way the 
Einstein field equations with singular sources. We therefore assume that the correct 
theory is general relativity in D = d + 1 space-time dimensions. As usual, any inter- 
mediate formulas will be interpreted by analytic continuation for a general complex 
spatial dimension d 6 C. In particular we shall analytically continue d down to the 
value of interest 3 and pose 

d = 3 + £. (89) 

The Einstein field equations in d spatial dimensions take the same form as pre- 
sented in Section [2] with the exception that the explicit expression of the gravita- 
tional source term A a P now depends on d. We find that only the last term in ( fTTT i 
acquires a dependence on d; namely the factor \ in g pa gen — jgaegpn should now 
read -pj. In addition the <f -dimensional gravitational constant is related to the usual 
three-dimensional Newton constant G by 



(90) 
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where £q is a characteristic length associated with dimensional regularization. 

In the post-Newtonian iteration performed in d dimensions we shall meet the 
analogue of the function F, which we denote by F^ (x) where x S W 1 . It turns out 
that in the vicinity of the singular points y^, the function FW admits an expansion 
richer than in (|82| |. and of the type 

F«(x) = f t^fttM + otf). (91) 



P=PO 9=90 



,-1 



The coefficients Afp.q (ni) depend on the dimension through £ = d — 3 and also on 
the scale £q. The powers of are now of the type p + qe where the two relative 
integers p,^eZ have values limited as indicated. Because reduces to F when 
e = we necessarily have the constraints 

( 92 ) 

9=90 A A 

To proceed with the iteration we need the Green function of the Laplace operator 
in d dimensions. Its explicit form is 

uf =Kr 2 A - d . (93) 

It satisfies Au^ = —4n8 A \ where 8^ = 8^ (x — y^) denotes the cZ-dimensional 
Dirac function. The constant K is given by 

K=^ T ^ 1 (94) 

where F is the usual Eulerian function. It reduces to one when d ^>3. Note that the 
volume Q c /_ i of the sphere with d — 1 dimensions is related to K by 

Q ^=w^- (95) 

With these results the Poisson integral of fW, constituting the <i-dimensional ana- 
logue of d86l) . reads 

In dimensional regularization the singular behavior of this integral is automatically 
taken care of by analytic continuation in d. Next we evaluate the integral at the 
singular point x' = y^. In contrast with Hadamard's regularization where the result 
was given by d87b . in dimensional regularization this is quite easy, as we are allowed 
to simply replace x' by y^ into the explicit integral form ( 1961 . So we simply have 
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d c/ x 

Z^2 



F^(x). 



(97) 



The main step of our strategy ll72l l73l will now consist of computing the differ- 
ence between the cZ-dimensional Poisson potential (l97l i and its 3-dimensional coun- 
terpart which is defined from Hadamard's regularization as (f8Tb . We shall then add 
this difference (in the limit e = d — 3 — > 0) to the result obtained by Hadamard reg- 
ularization in order to get the corresponding dimensional regularization result. This 
strategy is motivated by the fact that as already mentionned most of the terms do not 
present any problems and have already been correctly computed using Hadamard's 
regularization. Denoting the difference between the two regularizations by means of 
the script letter $1, we write 



®(P)a=P W (ja)-(P)a. 



(98) 



We shall only compute the first two terms of the Laurent expansion of @(P)a when 
£ -> 0, which will be of the form 2>{P)a = fl-i £ _1 + «o + This is the in- 

formation needed to determine the value of the ambiguity parameters. Notice that 
the difference &(P)a comes exclusively from the contribution of terms developing 
some poles <x l/e in the of-dimensional calculation. The ambiguity in Hadamard's 
regularization at 3PN order is reflected by the appearance of poles in d dimensions. 
The point is that in order to obtain the difference @(P)a we do not need the expres- 
sion of for an arbitrary source point x but only in the vicinity of the singular 
points y^. Thus this difference depends only on the singular coefficients of the local 
expansions of near the singularities. We find ||73l 



1 







(-+£ 




\9 





A H 



£=0 



,1+e 
AB 



■0(e). 



(99) 



We still use the bracket notation to denote the angular average but this time per- 
formed in d dimensions, i.e. 



A 



AO.. 



75 J pA^a) 

ii d-l A 



(100) 



The above differences for all the Poisson and interated Poisson integrals compos- 
ing the equations of motion (i.e. the accelerations of the point masses) are added to 
the corresponding results of the Hadamard regularization in the variant of it called 
the "pure Hadamard-Schwartz" regularization (see ll73l for more details). In this 
way we find that the equations of motion in dimensional regularization are com- 
posed of a pole part oc l/e which is purely 3PN, followed by a finite part when 
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£ — > 0, plus the neglected terms It has been shown (in the two-body case 

N = 2) that: 

1 . The pole part « 1 /e of the accelerations can be renormalized into some shifts of 
the "bare" world-lines by — > + | A , with % A containing the poles, so that 
the result expressed in terms of the "dressed" world-lines is finite when e — > 0; 

2. The renormalized acceleration is physically equivalent to the result of the Hada- 
mard regularization, in the sense that it differs from it only by the effects of shifts 
| A , if and only if the ambiguities in the Hadamard regularization are fully and 
uniquely determined (i.e. take specific values). 

These results l73l provide an unambiguous determination of the equations of motion 
of compact binaries up to the 3PN order. A related strategy with similar complete 
results has been applied to the problem of multipole moments and radiation field of 
point-mass binaries J74|. This finally completed the derivation of the general rela- 
tivistic prediction for compact binary inspiral up to 3PN order (and even to 3.5PN 
order). In later Sections we shall review some features of the 3.5PN gravitational- 
wave templates of inspiralling compact binaries. 

Why should the final results of the employed regularization scheme be unique, in 
agreement with our expectation that the problem is well-posed and should possess 
a unique physical answer? The results can be justified by invoking the "effacing 
principle" of general relativity 12011 — namely that the internal structure of the com- 
pact bodies does not influence the equations of motion and emitted radiation until 
a very high post-Newtonian order. Only the masses niA of the bodies should drive 
the motion and radiation, and not for instance their "compactness" Giua/ {c 2 cia). A 
model of point masses should therefore give the correct physical answer, which we 
expect to be also valid for black holes, provided that the regularization scheme is 
mathematically consistent. 



11 Energy and flux of compact binaries 

The equations of motion of compact binary sources, up to the highest known post- 
Newtonian order which is 3.5PN, will serve in the definition of the gravitational- 
wave templates for two purposes: 

1 . To compute the center-of-mass energy E appearing in the left-hand-side of the 
energy balance equation to be used for deducing the orbital phase, 

dF 

=- = -&', (101) 
at 

2. To order-reduce the accelerations coming from the time derivatives of the source 
multipole moments required to compute the gravitational-wave energy flux & in 
the right-hand-side of the balance equation. 
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We consider two compact objects moving under purely gravitational mutual in- 
teraction. In a first stage we assume that the bodies are non-spinning so the motion 
takes place in a fixed plane, say the x-y plane. The relative position X = yi — y2, 
velocity v = dx/df, and acceleration a = dv/df are given by 

x = rn, (102a) 
\ = rn + r(oX, (102b) 
a= (r-rco 2 )n+(r(b + 2r(o)X, (102c) 

The orbital frequency co is related in the usual way to the orbital phase (j) by ft) = 
(time derivatives are denoted with a dot). Here the vector X = z x n is perpendicular 
to the unit vector i along the z-direction orthogonal to the orbital plane, and to the 
binary's separation unit direction n = x/r. 

Through 3PN order, it is possible to model the binary's orbit as a quasi-circular 
orbit decaying by the effect of radiation reaction at the 2.5PN order. The restric- 
tion to quasi-circular orbits is both to simplify the presentational an d f° r physical 
reasons because the orbit of inspiralling compact binaries detectable by current de- 
tectors should be circular (see the discussion in Section [TJ. The radiation-reaction 
effect at 2.5PN order yieldQ 

f = -«^L v ?n +e (Ly (10 3a> 

*-f*»>* 1 +*(?). <"» 

where 7 is defined as the small [i.e. y = &{c~ 2 )\ post-Newtonian parameter 

Gm 

r= — - (104) 

rc 

Substituting these results into ( 1102| i, we obtain the expressions for the velocity and 
acceleration during the inspiral, 

y = r(a x- 6 Af9™ v ^ a+ ( }_] (105a) 

5 V r \c ' J 

■) 32 Gm en ( 1 \ 

a = -co 2 x-y J— vf l2 y+0[ ^jj . (105b) 

Notice that while r = (?(c~ 5 ), we have r = ^(c~ 10 ) which is of the order of the 
square of radiation-reaction effects and is thus zero with the present approximation. 



However the 3PN equations of motion are known in an arbitrary frame and for general orbits. 
13 Mass parameters are the total mass m = mi +ni2, the symmetric mass ratio v = mimi/m 1 
satisfying < v < 1/4, and for later use the mass difference ratio A = (mi — 1112) /m. 
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A central result of post-Newtonian calculations is the expression of the orbital 
frequency © in terms of the binary's separation r up to 3PN order. This result has 
been obtained independently by three groups. Two are working in harmonic co- 
ordinates: Blanchet & Faye ||69l l75l 1731 use a direct post-Newtonian iteration of 
the equations of motion, while Itoh & Futamase ll76l l77l |7§1 apply a variant of 
the surface-integral approach (a la Einstein-Infeld-Hoffmann 0) valid for compact 
bodies without the need of a self-field regularization. The group of Jaranowski & 
Schafer [79, 80 72] employs Arnowitt-Deser-Misner coordinates within the Hamil- 
tonian formalism of general relativity^ The 3PN orbital frequency in harmonic 
coordinates is 



Gm 
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(106) 



vMy 3 



Note the logarithm at 3PN order coming from a Hadamard self-field regularization 
scheme, and depending on a constant r'o defined by mlnr'o = milnr'i + m2lnr'2, 
where r^ = |x' — are arbitrary "constants" discussed in Section [9] We shall see 
that r'o disappears from final results — it can be qualified as a gauge constant. 

To obtain the 3PN energy we need to go back to the equations of motion for 
general non-circular orbits, and deduce the energy as the integral of the motion as- 
sociated with a Lagrangian formulation of (the conservative part of) these equations 
ll75ll . Once we have the energy for general orbits we can reduce it to quasi-circular 
orbits. We find 
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(107) 



A convenient post-Newtonian parameter x = G[c 2 ) 
defined from the orbital frequency as 



is now used in place of 7; it is 



Gmco 



2/3 



(108) 



The interest in this parameter stems from its invariant meaning in a large class of co- 
ordinate systems including the harmonic and ADM coordinates. By inverting ( 1106b 
we find at 3PN order 
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(109) 



This approach is extensively reviewed in the contribution of Gerhard Schafer in this volume. 
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This is substituted back into ( 1107b to get the 3PN energy in invariant form. We hap- 
pily observe that the logarithm and the gauge constant rg cancel out in the process 
and our final result is 
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The conserved energy E corresponds to the Newtonian, 1PN, 2PN and 3PN con- 
servative orders in the equations of motion; the damping part is associated with 
radiation reaction and arises at 2.5PN order. The radiation reaction at the dominant 
2.5PN level will correspond to the "Newtonian" gravitational-wave flux only. Hence 
the flying-color 3PN flux & we are looking for cannot be computed from the 3PN 
equations of motion alone. Instead we have to apply all the machinery of the post- 
Newtonian wave generation formalism described in Sections|4]-[7] The final result at 
3.5PN order is 
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(111) 



Here C = 0.577 • • • is the Euler constant. This result is fully consistent with black- 
hole perturbation theory: using it Sasaki & Tagoshi ll83~l [841 [851 obtain ( 111 lb in 
the small mass-ratio limit v — > 0. The generalization of (llllb to arbitrary eccentric 
(bound) orbits has also been worked out | 



12 Waveform of compact binaries 

We specify our conventions for the orbital phase and polarization vectors defining 
the polarization waveforms d52l in the case of a non-spinning compact binary mov- 
ing on a quasi-circular orbit. If the orbital plane is chosen to be the x-y plane as in 
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Section [TT] with the orbital phase measuring the direction of the unit separation 
vector n = x/r, then 

n = x cos + y sin0 , (112) 



where x and y are the unit directions along x and y. Following [88 89 1 we choose 
the polarization vector P to lie along the x-axis and the observer to be in the y-z 
plane in the direction 

N = y sin/ + z cos/, (113) 

where i is the orbit's inclination angle. With this choice P lies along the intersection 
of the orbital plane with the plane of the sky in the direction of the ascending node, 
i.e. that point at which the bodies cross the plane of the sky moving toward the 
observer. Hence the orbital phase <p is the angle between the ascending node and the 
direction of body 1. The rotating orthonormal triad (n, A,z) describing the motion 
of the binary and used in ( 1102b is related to the fixed polarization triad (N, P, Q) by 

n = P cos0 + (Q cos/ + N sin/) sin0 , (114a) 
X = -P sin</> + (Q cos/ + N sin/) cos0 , (114b) 
z = -Q sin/ + N cos/. (114c) 

The 3.5PN expression of the orbital phase as function of the orbital frequency 
or equivalently the x-parameter is obtained from the energy balance equation ( 1 1 1 b 
in which the binary's conservative center-of-mass energy E and total gravitational- 
wave flux J? have been obtained in (II 101 l- (ll 1 11 1. For circular orbits the orbital phase 
is computed from 

CO dE 
jFdco' 

Various methods (numerical or analytical) are possible for solving dl 15b given the 
expressions (II 10b and (II 1 lb . This yields different waveform families all valid at 
the same 3.5PN order, but which may differ when extrapolated beyond the normal 
domain of validity of the post-Newtonian expansion, i.e. in this case very near the 
coalescence. Such differences must be taken into account when comparing the post- 
Newtonian waveforms to numerical results l52l . 

It is convenient to perform a change of phase, from the actual orbital phase to 
the new phase variable 

. 2GMm ( co\ 
y/=</> 3 — In — , (116) 

where M is the binary's total mass monopole momen{3 an d &>o = 4^ ex P[i4 — C] 
is related to the constant uq = tq/c entering the tail integrals in < f6Qb . The logarith- 
mic term in \j/ corresponds physically to some spreading of the different frequency 
components of the wave along the line of sight from the source to the detector, and 



15 The mass monopole M differs from m = m\ +ni2 as it includes the contribution of the gravita- 
tional binding energy. At 1PN order it is given for circular orbits by M = m[l — ^y] + £?(c~ 4 ). 
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expresses the tail effect as a small delay in the arrival time of gravitational waves. 
This effect, although of formal 1.5PN order in (11 161 l, represents in fact a very small 
modulation of the orbital phase: compared to the dominant phase evolution whose 
order is that of the inverse of 2.5PN radiation reaction, this modulation is of order 
4PN and can thus be neglected with the present accuracy. 

The spherical harmonic modes of the polarization waveforms can now be ob- 
tained at 3PN order using the angular integration formula $55[ . We start from the 
expressions of the wave polarizations h+ and h x as functions of the inclination an- 
gle ; and of the phase l//. We use the known dependence of the spherical harmonics 
d54T > on the azimuthal angle. Denoting h = h + —ih x = h(i, y/) we find that the latter 
angular integration becomes 

//'" = (-i) m e - im v / d^' / di sin i h(i, y/) Y% (/, y/) , (117) 
Jo Jo 

exhibiting the azimuthal factor appropriate for each mode. Let us introduce 

a normalized mode coefficient H [m starting by definition with one at the Newtonian 
order for the dominant mode having (£,m) = (2,2). This means posing 



hlm= 2Gmvx /lftc 4. ^ 
Rc 2 V 5 

All the modes have been given in [54 1 up to the 3PN order. The dominant mode 
(2,2), which is primarily needed for the comparison between post-Newtonian cal- 
culations and numerical simulations, reads at 3PN order 

H22=1+ {- 1 -^ +5 i v ) x+2nx3/2 (119) 

2173 1069 2047 ,\ , / 107 34 „„. \ , /2 
-V 2 }x 2 +{ Tt-\ 7TV-241V x 5/2 



1512 216 1512 J \ 21 21 
/ 27027409 856 2 , 428 , s 428 
(-646800- I05 C+ r -T05 ln(lfo)+ K)5 i7r 



278185 41 , 

1 % 

33264 96 



20261 2 114635 ,\ 3 ^( 1 



2772 99792 / \c 



13 Spin-orbit contributions in the energy and flux 

To successfully detect the gravitational waves emitted by spinning, precessing bi- 
naries and to estimate the binary parameters, spin effects should be included in the 
templates. For maximally spinning compact bodies the spin-orbit coupling (linear in 
the spins) appears dominantly at the 1.5PN order, while the spin-spin one (which is 
quadratic) appears at 2PN order. The spin effect on the free motion of a test particle 
was obtained by Papapetrou 0641 in the form of a coupling to curvature. Seminal 
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later works by Barker & O'Connell |90l |9"T| obtained the leading order spin-orbit 
and spin-spin contributions in the post-Newtonian equations of motion. Based on 
these works, the spin-orbit and spin-spin terms were obtained in the radiation field 
Il92ll65ll93ll94ll , enabling the derivation of the orbital phase evolution (the crucial 
quantity that determines the templates). Finding the 1PN corrections to the leading 
spin-orbit coupling in both the (translational) equations of motion and radiation field 
was begun in 19511961 and completed in 119711981 . The result 1971 for the equations of 
motion was confirmed by an alternative derivation based on the ADM-Hamiltonian 
formalism l99l . 

In Section [8] we discussed a covariant formalism for spinning particles 115911601 
IBTl l62l l63l . We want now to find a convenient three-dimensional variable for the 
spin. Restricting ourselves to spin-orbit effects, i.e. neglecting ff(S 2 ), we can write 
the components of the spin tensor S A as 

Sf = - 1 — e ijk uj S A k , (120a) 

c^Wa 1 



*a = i^=e i]k 



A ™4 i A A oA 
U o k +U k — 0, 



C 



(120b) 



We have used the momentum-velocity relation d76l i and have taken into account 
the spin condition ([73}. A first possibility is to adopt as the basic spin variable the 
contravariant components of the spin covector Sf in (1120b , which are obtained by 
raising the index by means of the contravariant spatial metric, viz 

S A = (S A ) with S A = (Y j ) A Sj, (121) 

where is the inverse of the covariant spatial metric Yij = gij <•<?• satisfies / % ; = 
8j. The choice of spin variable (11211 i has been adopted in Il95ll96l . 

However to express final results (to be used in gravitational-wave templates) it is 
better to use a different set of spin variables characterized by having some conserved 
Euclidean lengths. Such spins with constant Euclidean magnitude will be denoted 
by S^. They can be computed in a straightforward way at a given post-Newtonian 
order in terms of the previous variables dl2U . For instance we find up to 1PN order, 



a A 



1 + ^ 



Sa-^(va-Sa)va + @[ - a ) • 



The (regularized) gravitational potential (U)a is defined by (ISTl ). The constant- 
magnitude spin variable obeys a spin precession equation which is necessarily 
of the form 



dS 
~df 



-G A xS A . (123) 



Indeed this equation implies that \S A \ = const. The precession angular frequency 
vector Q A for two-body systems has been computed for the leading spin-orbit and 
spin- spin contributions l65l and for the 1PN correction to the spin-orbit | 
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For two bodies we conveniently use the following combinations (introduced in |65l ) 
of the two spins: 



E c 



ni2 mi 



(124a) 
(124b) 



Furthermore, recalling the orthonormal triad {n,A,z} used in Section [TT] where z 
is the unit vector in the direction perpendicular to the orbital plane, we denote by 
SI = S c • z and E^ = E c -z the projections along that perpendicular direction. 

The spin-orbit terms have been computed at 1PN order both in the equations of 
motion and in the radiation field. In the equations of motion they correct the orbital 
frequency and invariant conserved energy with terms at orders 1.5PN and 2.5PN. 
In the presence of spins the energy gets modified to E = E mono 
monopole part has been obtained in (111 01 > and where the spin terms read 



£ sp in where the 



J spm 



2Gm 

fx 5 / 2 



■ is? 



-2AES 



10 



v)AEl 



(125) 



We recall that A = (mi — mij/m; see the footnote [T3l This expression is valid for 
quasi-circular orbits, and we neglect the spin-spin terms. Similarly the gravitational- 
wave flux will be modified at the same 1.5PN and 2.5PN orders. Posing & = 
^mono + ^spin where cFmono is given by (llllt . we find 



<3t ■ — 
■^spin 



spin 

32c 5 
5G 2 m : 

+ rV2 



v 2 x 5 



r 3/2 



-4S° 7 - -A Ei 
z 4 z 



272 



v )Si 



13 43 
-16 + 4 Vl ^ 



l-s . (126) 



Having in hand the spin contributions in E and we can deduce the evolution 
of the orbital phase from the energy balance equation dlOU . In absence of pre- 
cession of the orbital plane, e.g. for spins aligned or anti-aligned with the orbital 
angular momentum, the gravitational-wave phase will reduce to the "carrier" phase 
0Gw = 20 (keeping only the dominant harmonics), where is the orbital phase 
which is obtained by integrating the orbital frequency. However, in the general case 
of non-aligned spins, we must take into account the effect of precession of the or- 
bital plane induced by spin modulations. Then the gravitational-wave phase is given 
by $gw = 0gw + 50gw, where the precessional correction 50gw arises from the 
changing orientation of the orbital plane, and can be computed by standard meth- 
ods using numerical integration 01001 . Thus, the carrier phase 0gw constitutes the 
main theoretical output to be provided for the gravitational-wave templates, and can 
directly be obtained numerically from using the integration formula (11 15b . 
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